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Abstract 

The spin-polarization at the V(001) surface has been studied by using dif- 
ferent local (LSDA) and semilocal (GGA) approximations to the exchange- 
correlation potential of DFT within two ah initio methods: the all-electron 
TB-LMTO-ASA and the pseudopotential LCAO code SIESTA (Spanish Ini- 
tiative for Electronic Simulations with Thousands of Atoms). A comparative 
analysis is performed first for the bulk and then for a N-layer V(001) film 
(7<N<15). The LSDA approximation leads to a non magnetic V(001) sur- 
face with both theoretical models in agreement (disagreement) with magneto- 
optical Kerr (electron-capture spectroscopy) experiments. The GGA within 
the pseudopotential method needs thicker slabs than the LSDA to yield zero 
moment at the central layer, giving a high surface magnetization (1.70 Bohr 
magnetons), in contrast with the non magnetic solution obtained by means 

of the all-electron code. 

PACS numbers: 75.70.-i; 75.30.Pd; 73.20.-r 
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I. INTRODUCTION 



The understanding and prediction of the magnetic behavior of vanadium-based systems 
has motivated numerous experiments and calculations over the last fifteen years. Vana- 
dium is one of those paramagnetic metals! that can exhibit magnetism under certain condi- 
tions (loss of coordination, hybridization with a ferromagnet), due to its large paramagnetic 
susceptibility!'!. For instance, Akoh and Tasaki! have reported large localized magnetic mo- 
ments in hyperfine particles of V, and several experimental groups have demonstrated the 
existence of an induced magnetization at the V interface for V overlayers on Fe substratesIH! 
and Fe/V multilayers!~0. Although these trends are clear, in other aspects there has not 
been consensus. One controversial aspect concerns the short or long range induced spin- 
polarization in V. Two of us have reported on this problem in a recent wor where the 
reader can find a complete review. Another unclear aspect, that gave rise to an interesting 
discussion more than ten years ago, was the magnetic character of the V(001) surface. Rau 
et a0 through electron-capture spectroscopy concluded the existence of ferromagnetic order 
at the V(001) p(lxl) surface and on the V monolayer supported on Ag(OOl), whereas Fink et 
through magneto-optical Kerr measurements did not find magnetization in ultrathin epi- 
taxial films of V on Ag substrates. Ab initio full-potential linearized augmented plane wave 
(FLAPW) calculations within the density functional framework using the von Barth- Hedinlll 
local spin density approximation (LSDA) for the exchange and correlation (XC) potential 
found no surface magnetizationlll, which was since then admitted, and corroborated using 
other all-electron methods with LSDA (LMTO-ASA)El. 

Very recently, Bryk et aill have reported a theoretical work entitled "Magnetization of 
the V(001) surface in the generalized gradient approximation". They performed an ab initio 
calculation for a seven layer V(001) film using the plane-waves method with the ultrasoft 
(US) non-normconserving pseudopotential (PP) of VanderbiltEl, and the generalized gradi- 
ent approximation (GGA) of Perdewll (PW91) for the XC potential, obtaining a magnetic 
moment of 1.70//B at the unrelaxed surface plane ( 1.45/ie if relaxed). These authors be- 
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lieve that their prediction of a magnetic V(001) surface is more likely to be correct than 
the contrary result obtained earlierEl from local spin density approximation (LSDA), and 
consequently they request new measurements to test the prediction. Besides, although they 
correctly find the paramagnetic bulk ground state for V at the calculated equilibrium lattice 
constant, they also obtain, at the experimental lattice constant, the ferromagnetic phase, 
with 0.15/iB P er atom, lying only 0.7 meV below the paramagnetic one. In the opinion 
of the authors, this energy is small compared with the difference between the calculated 
magnetic and non magnetic surface energies, concluding that the GGA overestimates the 
magnetization of the interior layers, but does not alter their prediction of magnetic V(001) 
surface. An even more recent calculation of the magnetic structure of V(001) surface using 
the all- electron FLAPW method and the PW91 form of GGA§, by Bihlmayer, Asada and 
Bliigel@, conclude that "in very thin V films a surface magnetic moment can be stabilized, 
while for thicker and relaxed films no surface magnetism can be found" . 

The suspicion that the PW91 form of GGA might incorrectly tip metals with large LSD 
spin-susceptibility enhancements, like Pd and V, into the ferromagnetic state, was posed in 
the first paper testing the PW91 functional^. Subsequent calculations using all-electroniHl 
and/or pseudopotentialsiiHil methods, coincide pointing out the trend of GGA's to enhance 
magnetism of magnetic materials and susceptibilities of non magnetic materials when com- 
pared with results from analogous LSDA calculations. Singh and AshkenaziH using the 
all-electron FLAPW method confirmed that bulk V and Pd are correctly predicted as para- 
magnetic by the PW91 GGA approximation, but conclude also that GGA's do not have 
greater precision than the LSDA for studying transition metals (TM) and specially for 
magnetic materials. In Ref. [20] it is also noted that there are considerable differences 
between the predictions of the GGA's from Langreth-Mehl-Hu (LMH)0 and PW91, which 
gives room for more consistent improvements on the LSDA. Ozolins and Korlingil, using 
the full-potential linear muffin-tin-orbital (FP-LMTO) methocSi found that the PW91 
GGA predicts more accurately than the LSDA the structural properties of non magnetic 
3d, 4d and 5c? TM. Also in Ref. [21] it is suggested that a full-potential treatment, instead 



of the atomic-sphere approximation (ASA), is important for GGA calculations, because the 
inaccuracies introduced by the ASA are of the same magnitude as the gradient corrections. 
Apparently, an ASA GGA calculation average out the angular parts of the GGA and thus is 
missing what could be a portion of the GGA contribution. Norm conserving and ultrasoft 
pseudopotentials has been compared for first-row and transition elements by Kresse and 



transferability properties. The extreme softening of non-normconserving Vanderbilt pseu- 
dopotentials without loss of accuracy compared to very hard and accurate pseudopotentials 
(working at much more higher cutoff energy) has been demonstrated recently by Furthmuller 
and coworkers!!. On the other hand, the use of pseudopotentials, particularly the US-PP of 



Whereas Holzwarth et al" conclude that the structural properties of bcc V are represented 
equally well by the PAW, LAPW and PP methods, Kresse and Joubert@ point out that 
for TM at the left side of the Periodic Table it is very desirable to include semicore states 
as valence states in the pseudopotentials. The LSDA (GGA) magnetic moments of bcc Fe 
calculated in Ref. [23] are slightly (considerably) larger within the PP method than within 
the AE methods FLAPW and PAW. On the other hand, in both PP and AE calculations, 
the GGA magnetic moments of Fe, Co and Ni are systematically larger than the LSDA ones. 
The same trends are obtained in Ref. [24] by comparing calculations using the PP-LCAO 
code SIESTA§ with available AE results for the magnetic moment of bulk bcc Fe, the (001) 
Fe surface, and small clusters of Fe on the (001) Ag surface. 

The work of Bryk et a0 opens again clearly the discussion about the magnetic character 
of V(001) and besides, from the theoretical point of view, it offers a new benchmark to test 
different ab initio methods and different XC functionals. Therefore, the need of further cal- 
culations is evident. In this work, calculations have been done using two different ab initio 
methods based on the density functional theory (DFT)il, namely, i) AE tight-binding linear 




from the point of 



view o numerical performance achieving good convergence and 
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muffin tin orbitals (TB-LMTO) method, in the scalar relativistic version and the ASA§§E3, 
and ii) the SIESTA LCAO methodic, using the Troullier-Martins pseudopotentialsll. For 
the XC functional, we have used within TB-LMTO four different approximations: two LSDA 
versions -the Hedin-Barth0 and the Vosko-Wilk-Nusairii-, and two GGAs -the PW91@ and 
the Langreth-Mehl-Hu (LMH)H. The XC functionals used within the pseudopotential code 
SIESTA0 are the Perdew-Zunger (PZ) parameterization for the LSDA0 and the Perdew- 
Burke-Ernzerhof (PBE)@ form of GGA. The soft ionic pseudopotential used in SIESTA is 
generated according to the procedure of Troullier and Martinsi! from the atomic configura- 
tions [Ar]3d 3 4s 2 for V with core radii for the the s, p and d components of 2.35, 2.70 and 
2.35 a.u. respectively. The partial core-correction for non-linear XC§ has been included. A 
careful study of the optimum core- correct ion radius leads to a value of 0.8 a.u. 

Let us first discuss the bulk bcc case. The GGA's equilibrium lattice constants are, 
within the TB-LMTO, (in a.u.) 5.68 (5.79) from the LMH (PW91) calculations, and 5.71 
for the SIESTA-PBE calculation. The LSDA's equilibrium lattice constant is 5.61 a.u. for 
the different LSDA approximations in both TB-LMTO and SIESTA calculations. With 
both AE and PP methods, and for all the XC-approximations used we obtain a paramag- 
netic ground state at the experimental lattice constant (5.73 a.u.), in agreement with the 
experimental findings^!. This contrasts with the GGA result of Bryk et al3 where bulk V at 
the experimental lattice constant is ferromagnetic (F) with 0.15/^s per atom, although they 
find the paramagnetic state at the equilibrium lattice constant in their calculations. It is a 
general trend that when expanding the lattice, the electronic localization increases and the 
kinetic energy of the system decreases so that spin-polarization is favored. It is therefore ex- 
pected to find a magnetic phase transition in V bulk in this context. Moruzzi and Marcus^ 
studied ten years ago this magnetic transition in V bulk through total-energy band calcu- 
lations using the LSDA in the augmented-spherical-wave (ASW) method. They found that 
the paramagnetic solution is the ground state of V bulk for the experimental lattice constant 
and for lattice expansions of less than 12%. For this expansion they find a transition to the 
antiferromagnetic (AF) solution (always more stable than a low-spin F solution). Earlier 
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spin-polarized augmented plane wave (APW) calculations for bcc vanadium! found also that 
a non magnetic to ferromagnetic transition occurs abruptly for a lattice constant about 25% 
larger than the equilibrium value. We have also tested volume expansions and within all 
XC-functionals considered in our TB-LMTO and SIESTA calculations, the paramagnetic 
ground state persists at least up to an expansion of 10%. 

In the case of the seven layer V(001) film at the bulk experimental lattice constant, the 
LSDA gives no surface magnetization in both our AE and PP calculations, in agreement 
with previous first principles LSDA result JiE. Concerning the GGA, our all electron results 
are reported in figure 1 for both the PW91 and LMH versions of GGA. Here a surface 
magnetization is obtained, although it is much lower than the value reported by Bryk et a0 
using pseudopotentials. They obtain 1.70/z.b at the surface plane in contrast with our values 
of 0.66/iB for the PW and 0.25/is for the LMH cases (Fig. 1). Notice also the large difference 
obtained with different versions of GGA, as already observed in Ref . [20] , and that for volume 
expansions of the order of 2% that difference become negligeable (Fig. lb). We see then, 
that the AE surface magnetization is still far from the PP value of 1.1 Hb obtained by Bryk 
et a£3. This system constitutes a good test for an ab initio pseudopotential calculation, 
in order, to see whether or not the pseudopotential code give rise to an overestimation of 
the V magnetic moment compared with a typical AE method. For this purpose we have 
repeated the same seven-layer calculations with SIESTA and the PP of Troullier-Martins. 
The LSDA again leads to the paramagnetic solution whereas the GGA leads to a magnetic 
moment of 1.77 '//g which is similar to the value reported by Bryk et alii. This results 
indicate that the GGA within the pseudopotential calculations enhances even more than 
within AE methods the V (001) magnetic moment. We note, however, a difference between 
our PP-GGA results and those of Bryk et for the seven layer V film. At the experimental 
lattice constant, they achieve in the central layer the convergence of the magnetic moment to 
their bulk value (0.15//b) and, consequently, they use the surface layer to discuss the V(001) 
surface. However, our results for the central layer do not converge to our paramagnetic 
bulk with just seven layers. As, we increase the number of layers in the film, the calculated 



magnetic moment at the center decreases. For 15 layers the local moment at the center 
is zero as in the bulk and we can discuss then the V(001) surface magnetism in terms of 
the surface layer. The interesting result is that when we consider more than seven layers 
in the slab the V(001) surface is non magnetic with the all-electron TB-LMTO method, 
independently of the approximation used for the XC potential (LDA or GGA), whereas 
with the pseudopotential code SIESTA, even for calculations with 15 layers of V the GGA, 
still produces a high magnetic moment at the surface (see table 1), comparable with that 
obtained by Bryk et a 10 (with LDA we obtain the non magnetic surface). 

In resume our results indicate the well known tendencyilil of GGA to produce larger 
magnetic moments than the LSDA, but also that the convergence of GGA to the bulk value 
at the central layer of a slab is slower than for LSDA. Taking care of these convergence 
problems, the V(001) surface is non magnetic with both the LSDA and GGA within the 
all-electron TB-LMTO method, whereas the pseudopotential code with GGA still produces 
a high magnetic moment. It is interesting to note that our pseudopotentials and basis 
in SIESTA are different from those of Bryk et al0. On the other hand, in view of the 
controversial theoretical predictions, new experiments concerning the magnetic character of 
the (001) surface of vanadium are needed to confirm or to correct the earlier findings of Rau 
and coworkers^. 
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FIGURES 

FIG. 1. Magnetic moments profile for the 7-layer V(001) film obtained with TB-LMTO method 
and two GGA versions (LMH0 and PW9lll!). Upper panel corresponds to the experimental lattice 
constant whereas lower panel corresponds to a lattice expansion of 2%. 
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TABLES 

TABLE I. Local magnetic moment obtained at the surface of a 7-layer and 15-layer V(001) 
film with the different versions of the GGA in the all-electron TB-LMTO and pseudopotential 
SIESTA methods for the experimental bulk lattice constant of vanadium. For slabs thicker that 7 
layers the TB-LMTO yields non magnetic ordering at the V(001) surface. 





TB-LMTO (LMH) 


TB-LMTO(PW91) 


SIESTA(PBE) 


7 layers 


0.25 


0.66 


1.77 


15 layers 


0.00 


0.00 


1.70 



12 



Magnetic Moment (|Li B ) 

o p 
OOI O Ul 





1 I I 




1 
1 


1 


□ 

II 


1 1 1 1 1 1 1 


! -S - 

:fl si 



